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Abstract. Simultaneous Buda-Lund hydro model fits are presented to iden- 
tified particle spectra and two-particle Bose-Einstein correlations as measured 
by the STAR collaboration in ^ = 200 GeV p+p collisions at RHIC. Pre- 
liminary results are compared to similar results in Au+Au collisions at RHIC, 
hadron-|-p reactions as well as Pb-|-Pb collisions at CERN SPS. 

Keywords: identified particle spectra, Quark-Gluon Plasma, Bose-Einstein 
(HBT) correlations, hydrodynamical models, soft proton-proton collisions 
PACS: 24.10.Nz, 25.75.-q, 25.75.Nq, 25.75.Gz 



X 



1. Introduction 

The Buda-Lund hydrodynamic model [^ |2] was shown to describe in a statistically 
acceptable manner the single particle spectra and the two-particle Bose-Einstein 
correlations in hadron-|-proton collisions at CERN SPS [|3]. The model, with differ- 
ent source parameters, also described well the identified single particle spectra and 
the two-particle Bose-Einstein correlations in Pb-|-Pb collisions at CERN SPS [^. 
Similarly, the Buda-Lund hydro model described well the single particle spectra and 
the two-particle Bose-Einstein correlation functions in Au-|-Au collisions at RHIC [ 
ISj, both at ^JSNN = 130 and 200 GeV colliding energies. The model was recently 
reviewed in ref. [E] and it was extended to describe elliptic flow in non-central heavy 
ion collisions at RHIC energies in refs. [[3 El- 

We have observed in these works, that the RHIC Au+Au data are well described 
by the model with a relatively high central temperature of 214 ± 7 MeV in 0- 
10 % central Au-hAu collisions at ^/sWn = 130 GeV and 200 ± 9 MeV central 
temperature in the less central, 0-30 % Au-f Au collisions at ^snn = 200 GeV. 
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The same analysis suggests, that the surface (or with other words, the average) 
temperature of the hadronic final state is remarkably low, around To/2, nearly 100 
MeV in these reactions. 

Lattice QCD calculations estimated the critical temperature of the hadron - 
quark gluon plasma phase transition near /x^ = to be Tc = 172 ± 3 MeV [^, 
which was recently calculated more precisely, using the physical quark masses, to 
be Tc = 162 ± 2 MeV [HHI, even lower than thought before. Thus the Buda-Lund 
data analysis of central Au+Au spectra and correlation radii implies the existence 
of a superheated hadron gas in the very middle of the Au+Au collisions at RHIC, 
providing an indirect indication for the transition to deconfined matter [|S]. How- 
ever, at the lower CERN SPS energy, the central temperature was found, in similar 
frameworks, to be within errors the smaller 139 ± 6 MeV, both in hadron+proton 
collisions [ Oj and in Pb+Pb collisions [ El , in line with the arguments of Landau 
that the freeze-out temperature of a hadron gas should be approximately the ir^ 
mass, the mass of the lightest neutral quanta that can be excited in such a system. 

In this work we investigate two questions: 

a) Can a hydrodynamical model describe successfully the final state (identified 
particle spectra and Bose-Einstein correlations) in p+p collisions at RHIC? 

b) What are the best hydro model parameters and how they correspond to the 
values for hadron-|-p reactions at CERN SPS and Au+Au collisions at RHIC ? 

To answer these questions, we utilize as a tool the Buda-Lund hydrodynamical 
model. This model has been successfully applied to describe the double-differential 
invariant momentum distribution of charged particles in hadron+proton [ T and 
Pb+Pb collisions [ 4 at CERN SPS fixed target experiments, as well in Au+Au 
collisions at the RHIC collider energies of ^/SNN — 130 GeV and 200 GeV [[5]. Very 
few, if any, other models pass all of these tests (as far as we know). 

2. The Buda-Lund hydro model 

The Buda-Lund hydro model was formulated in refs. [HE]) based on the following 
four principles: 

- The particle emitting source consist of two, physically very different regions, 
the core and the halo. The core is a fireball like a star, the halo consists of the decay 
products of long-lived resonances, with relatively large length-scales, reminiscent to 
the solar wind surrounding the star. The large length-scales in the halo are assumed 
to be larger than the maximal length-scale resolvable by particle interferometry. 

- The core is described by a locally thermalized, three dimensionally expanding 
particle emitting source. For central collisions, we utilize axial symmetry, for semi- 
central collisions, ellipsoidal symmetry to constrain the structure of the flow velocity, 
the fugacity and the local temperature distributions. 

- We parameterize the hadronic final state in such a way, that it corresponds 
to know solutions of (relativistic or non-relativistic) hydrodynamical solutions in 
certain domains of the parameter space. In these limiting cases, even the time 
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evolution of the local densities and the velocity field can be followed. 

- We constrain the parameterization in such a way, that the observables can 
be calculated analytically from the model. This is useful, as simple scaling laws 
are observed experimentally in the transverse momentum spectra as well as in the 
Bose-Einstein (HBT) correlation functions. The Buda-Lund hydro model observes 
these scaling behaviors in certain domains of the parameter space. 

2.1. Hydro in soft hadron+hadron collisions 



BudaLund v1 .5 fits - p+p data at 200 GeV 



BudaLund v1 .5 fits - p+p prel. at RHIC 
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Fig. 1. The Buda-Lund hydro model is fitted here simultaneously to STAR single 
particle SDectra[ llll I12j and STAR preliminary two-pion correlation radius parame- 
ters [El in P+P collisions at -^/s = 200 GeV, with small amount of transverse flow, 
but with a large central temperature Tq > Tc- The best fit parameters are summa- 
rized in Table I. The upper plot on the right side indicates the predicted pseudo- 
rapidity distribution of charged particles these collisions: BRAHMS and PHOBOS 
data can be used later, to restrict the source parameter Ary. 



Recently, the STAR collaboration reported the identified particle spectra and 
the two-particle Bose-Einstein correlation radii in proton+proton collisions at the 
maximal RHIC colliding energy of ^/snn = 200 GeV. At the Quark Matter 2004 
conference, the transverse mass dependence of these radius parameters was pre- 
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sented as one of the new HBT puzzles at RHIC [ 15 . However, the transverse 
mass dependence of the HBT radii in p+p coUisions at RHIC is very similar to the 
transverse mass dependence of these radii in hadron+proton reactions at the lower 
CERN SPS energies [EIj which were well described in terms of the Buda-Lund 
hydrodynamical model, together with the rapidity and transverse mass dependent 
double differential spectra of charged particles [ 3 . In that case, the transverse 
mass dependence of these HBT radii was not generated by the transverse flow, but 
by the transverse temperature inhomogeneities of hadron+proton collisions. 

Perhaps it is surprising, that hydrodynamical descriptions work in the soft 
domain of high energy hadron+hadron collisions. However, Fermi and Landau 
initiated the statistical and hydrodynamical description of particle production for 
nucleon-nucleon as well as for nucleus-nucleus collisions. After accepting the hydro- 
dynamical language, it is not surprising that the transverse temperature gradients 
play major role in hadron+hadron collisions - these systems have a rather small, 
typically 1 fm transverse radius, and the hot hadronic matter inside is surrounded 
with a vacuum that has zero temperature, so the temperature decreases on a small 
transverse radial scales, as compared to high energy heavy ion reactions, where the 
temperature decreases on a factor of 5-10 times larger scales. 

2.2. Buda-Lund fit results 

In this section we show fit results to STAR final identified particle spectra [El ^] 
and two-pion Bose-Einstein correlation radii [ E] in p+p collisions at -^/i = 200 
GeV. First we investigate the question, is it possible to fit these data with the same 
central freeze-out temperature and transverse flow parameter, as were obtained [E] 
in hadron+p colhsions at CERN SPS? 

The answer is no, as also show graphically in the presentation version of this 
work [Ell however, this bad fit cannot be shown here due to space limitations. 
But the conclusion is clear: When fixing the freeze-out temperature in the center 
of the fireball to Tq — 140 MeV and the transverse fiow to (ut) — 0.2, as in ref. [ 
E] , it is not possible to tune the other parameters of the Buda-Lund hydro model 
to describe the STAR p+p data at RHIC. 

However, the Buda-Lund hydro model can successfully describe all these ob- 
servables in p+p collisions, when the central temperature and the transverse fiow 
parameters are released, as shown in Fig. 1. From this comparison one learns, that 
the central temperature in p+p collisions at RHIC is significantly larger, than in 
h+p collisions at CERN SPS, although the transverse flow parameter is in both 
cases within errors compatible with zero. The fit parameters are summarized in 
Table 1. 

3. Conclusions 

Let us quote Landau, who wrote [ I17| in 1956: "Experiment shows that in collisions 
of very fast particles a large number of new particles are formed in multi-prong 
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parameter 


p+p 200 GeV 


Au+Au 


200 GeV 


Au+Au 


130 GeV 


To [MeV] 
Te [MeV] 
MB [MeV] 


289 
90 

8 


±8 
± 42 
± 76 


200 

127 

61 


± 9 
± 13 
± 40 


214 

102 

77 


± 7 
± 11 
±38 


Rg [fm] 
Rs [fm] 

K) 


1.2 

1.1 

0.04 


±0.3 
±0.2 
±0.26 


13.2 

11.6 

1.5 


± 1.3 
± 1.0 
± 0.1 


28.0 
8.6 
1.0 


± 5.5 
± 0.4 
±0.1 


To [fm/c] 
Ar [fm/c] 
A?/ 


1.1 
0.1 
3.0 


±0.1 
±0.5 
fixed 


5.7 
1.9 
3.1 


± 0.2 
± 0.5 
± 0.1 


6.0 
0.3 

2.4 


±0.2 
± 1.2 
±0.1 


XVNDF 


89.7 


/69 


132 


/208 


158.2 


/ 180 



Table 1. Buda-Lund hydro model vl.5 source parameters, corresponding to Fig. 
1. The errors and the fit parameters are preliminary, as some of the fitted points 
are not yet final and the statistical and systematic errors were added in quadrature 
in these fits. 



stars. The energy of the particles which produce such stars is of the order of 10"'^^ 
eV or more. A characteristic feature is that such collisions occur not only between 
a nucleon and a nucleus, but also between two nucleons." 

The Buda-Lund hydro model describes in a statistically acceptable manner the 
soft identified particle spectra and two-particle Bose-Einstein correlation data of 
STAR in p+p collisions at -^s = 200 GeV. The temperature distribution has a 
maximum which is significantly above the critical Tc = 162 ± 2 MeV value [llOj. 
but falls rapidly, within about 1.2 fm transverse distance, towards the vanishing 
temperature of the vacuum. 

Our results are in contrast to the recent statement by Shuryak [tlSj, who claimed 
that a hydrodynamical approach cannot describe the data in p+p collisions at the 
RHIC energy domain. It seems to us that Landau's insight was more precise in 
this particular case. Even more surprising, the flow velocity profile, within the 
errors of reconstruction, corresponds to the Hwa-Bjorken famous 1+1 dimensional 
relativistic hydrodynamical solution [ I19L I20| . In the transversally more extended 
Au+Au collisions we find evidence for a fully developed, three dimensional Hubble 
flow [0. For similar resuhs, see [ EH EH E31 E1I23 • 
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